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The use of logarithmic frequency spacing brings several advantages to the reduction of 
ionograms to electron density profiles. Among them is the fact that, when computing factors 
for the analysis, one need not determine the group refractive index. Formulas involving 
only the phase refractive index are presented; for the ordinary component one exact, and one 
approximate formula are given, while for the extraordinary component there is an approxi- 
mate formula valid over a wide range of geomagnetic latitudes. There is a brief discussion 
of quasi-longitudinal approximations to the extraordinary phase refractive index. 



1. Introduction 

In his important paper on obtaining electron den- 
sity profiles, Budden [l] 3 used for illustration a 
linear sampling of the frequency scale of the iono- 
gram. He also mentioned the possibility of other 
samplings, including logarithmic. At about the same 
time, King [2] presented an illustration based on 
logarithmic sampling, and later [3] showed that this 
spacing simplifies inclusion of the earth's magnetic 
field in the analysis, by the use of an accurate 
approximation. 

The methods of real height analysis using sum- 
mation of spaced ordinates due to Kelso [4], Shinn 
(described by Thomas [5]), and Schmerling [6], are 
much easier to apply if the ionograms have a loga- 
rithmic frequency scale, for then overlays can be 
used. This derives from the fact that the properties 
of the propagation equation depend mainly on Jn\J, 
the ratio of the plasma frequency to the exploring 
frequency, and to a much less extent on the absolute 
frequencies (determined byf H /J, where f H is the gyro- 
frequency. As is usual in treating this subject, the 
effects of electron collisions are neglected.) 

The purpose of this note is to inquire more closely 
into the advantages of logarithmic spacing and to 
present an approximation using it for analysis of the 
extraordinary ionogram trace. 

2. General Equations 

The basic equation relating the virtual height, h', 
and the real height, h n at the reflection point (as- 
suming geometrical optics) is, [2] 
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where y! is the group refractive index. This can be 
rewritten 



'=f*V- 



h'= 






(2) 



where <E> is any single valued function of the electron 
density. The integration can now be divided into a 
series of steps within which dh/d& is varying suffi- 
ciently slowly to be taken out of the integral sign. 

/4=S(~) f n'd* (3) 

m \A0/ m J A<t>m 

The integral n'd$) which we shall denote by F m 

defines constant factors used in the analysis of the 
ionograms to obtain h. Once the table of factors is 
prepared, the reduction to real heights is simply a 
matter of solving a set of simultaneous equations in 
(AA/A*) m . 

If now $ is identified with the logarithm of j Nj a 
simplification results. (The natural logarithm, In, 
will be discussed here to simplify the presentation, 
although common logarithms may be used in 
practice.) 

For, at fixed f N , 






i/a 



=M+ [anr/l/ 



(4) 



where \x is the phase refractive index; it is a function 
of both/ and /y. 
Using the theorem 



Lain/I^ 



"bto/l* + UL/_l 



(5) 
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relating the partial differential at fixed f N /f to partial 
differentials at fixed f N and at fixed/, we get 

The factors used in the ionogram analysis then 
become 



r M , ^= f 



1q /at 



lx f dh\j N 



^In/WA^+Dy^Aln^, 



(7) 



where the "bars" denote mean values of the quan- 
tities over the interval. 

All terms of eq (7) can be evaluated from a table 
of ^(ln/jln/jv), so that one can avoid the tedious 
calculation of the group refractive index. In a 
recent paper, Titheridge [7] has given an expression 
with the same advantage. 

3. Ordinary Component 

While one would normally use the complete eq 
(7), it is worth noting that, for routine analysis of 
the ordinary component, the last term can be 
dropped [3]. Experience shows that the heights 
deduced from the ionograms are not seriously affected 

The assumption here is 



[^1 



_dln/J /jv// 
4. Extraordinary Component 



(8) 



While eq (7) applies formally to either the ordi- 
nary or the extraordinary component its use with 
the latter is made difficult by the rapid change of 
the last term near the gyrofrequency. By making 
the transformation 



?=f(l-y), 



(9) 



where 2/ =/#//, one can adopt an assumption similar 
to eq (8), 

0, (10) 



r-^-i - 



and so simplify the computations. 

The approximation of eq (10) will be considered 
in the next section where it will be shown to hold 
extremely well down to geomagnetic latitudes as 
low as 30° (where the propagation angle, 0, between 
the earth's magnetic field and the vertical is as 
large as 50°). 

When eq (7) is put in terms of £ and eq (10) is 
applied, one obtains 

F xm = I fx'd In J N =jlA In f N ~j^[A»] v (11) 
jAin/jv am j 

where the subscript x in F Xi 
nary component. 



denotes the extraordi- 



As constant £ implies constant/, and 

dnf l-y/2 
dlnf 1-y ' 



(12) 



the factors for analysis of the extraordinary compo- 
nent become 



F xm =fx Aln f N 



i-y/2 

1-2/ 



[Am]/- 



(13) 



5. Phase Refractive Index for the Extra- 
ordinary Mode 

In order to justify use of the approximation (10), 
a few remarks on quasi-longitudinal approximations 
to the phase refractive index of the extraordinary 
mode are appropriate. 

Common approximations are [8] # 



(14) 





1 

1-M 2 


f 2 
In 


cos 6) 


id [9], 










1 
1-M 2 * 


P 
= /l(l-|f) = 

JN 


-m- 


The exact 


expression is 





(15) 



^=£i(l-2/cos<?cot|), (16) 



where 



tan ip— 



2(f-f N 2 ) cosfl 



SSh 



sin 2 6 



(17) 



Equations (14), (15), and (16) can be compared as 
follows: For the extraordinary trace, the least value 
of tan ^ occurs at the reflection point, where it is 
2 cos 0/sin 2 6, so that the greatest value of cot \///2 
is 1/cos 0. 

That is to say, 



l<cot^<-i- 

2 cos 6 



(18) 



Therefore the exact value for 1/1-m 2 (16) always 
lies between those given by the two approximate 
expressions (14) and (15), except at the reflection 
point, where (15) is exact. 

Now if we let A represent the ratio of (16) to (14), 
and B the ratio of (16) to (15), i.e., 



and 



A= 



B= 



1— y cos 6 cot ; 
1 — 2/ cos 



1— ycosd cot ; 
1-2/ 



(19) 



(20) 
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the following table compares the ratios for y—\ and 
various values of (J N [f) 2 , at different geomagnetic 
latitudes. 



Latitude 


df> 


0.5 


0.4 


0.3 


0.2 


0.1 


0.0 


14° 


A 
B 


0.644 
1.000 


0.714 
1.109 


0.762 
1.184 


0.797 
1.238 


0. 823* 
1.278 


0. 844* 
1.310 


23§° 


A 
B 


0.743 
1.000 


0.792 
1.065 


0.825 
1.110 


0.850 
1.143 


0. 868* 
1.168 


0. 883* 
1.187 


30° 


A 
B 


0.804 
1.000 


0.840 
1.045 


0.864 
1.075 


0.883 
1.098 


0. 897* 
1.116 


0. 908* 
1.129 


45° 


A 
B 


0.904 
1.000 


0.921 
1.018 


0.933 
1.031 


0.941 
1.041 


0.948 
1.048 


0. 954* 
1.054 


60° 


A 
B 


0.962 
1.000 


0.969 
1.007 


0.973 
LOU 


0.977 
1.015 


0.979 
1.018 


0. 982* 
1.020 



Only in those positions of the table marked by an 
asterisk is A closer to unity than B. Clearly, then, 
the approximation (15) is the better, especially near 
the reflection point [{j N lj) 2 = 1— ?/] where the need for 
accuracy is greatest. Accepting eq (15), eq (10) 
follows. 

Writing the exact expression, using eqa (16), (20), 
and (9), as 



1 -? 7? 



(21) 



it can be shown that, because of the slow rate of 
change of B with (J N /J) 2 (and hence with ?lf N 2 ) : eq 
(10) is applicable to ionogram analysis over a wider 
range of latitudes than eq (15). Ryclbeck [9] used a 
relation equivalent to eq (15) in a method of ionogram 
analysis for 0<2O° (latitudes greater than 54°), 
where it certainly holds within the accuracies to 
which the ionograms can be read. Use of eq (10) 
extends the range of application down to geomagnetic 
latitude 30°. 

6. A Numerical Example 

The following tables illustrate the computation of 
the factors F xm for analyzing the extraordinary 
ionogram trace. 

The first table gives n as a function of log 10 f N and 
logio £• This was calculated using eq (16). 



Table of /xX 10* 


f H = 1.4753 Mc/s 


0=21°53' 


f 


log* 


lOg fN 


Mc/s 


0.32 


0.28 


0.24 


0.20 


0.16 


0.12 


0.08 


0.04 


0.00 


2.95 
2.78 
2.63 


0.32 
0.28 
0.24 
0.20 
0.16 
0.12 
0.08 
0.04 
0.00 


0000 


4209 
0000 


5658 
4210 
0000 


6610 
5660 
4211 
0000 


7302 
6613 
5662 
4212 
0000 


7827 
7305 
6616 
5664 
4214 
0000 


8236 
7830 
7308 
6619 
5667 
4215 
0000 


8560 
8239 
7833 
7311 
6622 
5669 
4217 
0000 


8821 
8562 
8242 


2.49 






7836 


2 36 








7314 


2.25 










6625 


2 15 












5672 


2 06 














4219 


1 98 
















0000 























Lines parallel to the diagonal in the table define 
values of \x for constant values of the ratio /#/£. 
Differences between the values along such a line, 
therefore, give [dfx/d log %\ fN /$ ; and as this is very small 
compared with [d/x/d log £\ fiv and [d/z/d log J N ]^ the 
assumption made in section 4 (eq (10)) is quite 
satisfactory. 

From this table one can obtain the components 
mA \ogf N and M\(l — y/2)/(l— y)] [Am]/, of the factors 
Fxm (eq (13)). Here, M=0.4343 is the conversion 
factor to change from natural logarithms to common 
logarithms. Thus for log £=0.24, we obtain the 
factor for the interval A log /at (0.16 — 0.12) by adding 

Mlog^ (0.0240) and -M^^ [An] f (0.0680) giving 

a value of 0.0926. 

The table of F xm Xl0* is then, 





Interval of log fy 


log* 


0.28 

to 

0.32 


0.24 
to 
0.28 


0.20 

to 

0.24 


0.16 

to 
0.20 


0.12 
to 
0.16 


0.08 
to 
0.12 


0.04 

to 

0.08 


0.00 
to 
0.04 


-0.04 

to 
0.00 


-0.08 

to 
-0.04 


-0.12 

to 
-0.08 


0.32 
.28 
.24 


2855 


1142 
2976 


866 
1184 
3117 


730 
894 
1233 
3280 


645 
749 
926 
1290 
3470 


587 
660 
772 
964 
1356 
3692 


574 
599 
677 
799 
1006 
1434 
3954 


518 
556 
613 
697 
830 
1058 
1525 
4262 


494 
525 
566 
628 
721 
867 
1117 
1631 
4625 


476 
500 
534 
579 
646 
749 
910 
1187 
1757 


462 
481 
506 


.20 






544 


.16 








593 


.12 










668 


.08 












780 


.04 














960 


.00 
















1269 























This example is taken from part of a table prepared 
for testing methods of ionogram analysis which use 
both the ordinary and extraordinary traces to give 
information on the unobserved parts of the iono- 
sphere [10], [11], [12]. 

7. Conclusions 

The use of logarithmic frequency spacing allows 
easy computation of factors for the real height 
analysis of ionograms, using tables of phase refractive 
index; the group refractive index need not be calcu- 
lated. 

For the ordinary wave component, there is an 
exact formula for the factors eq (7), and an approxi- 
mate formula good enough for most work. The 
formula for the extraordinary wave component eq 
(13) contains an approximation, but it can be used 
with negligible error for analysis of ionograms from 
any but the lowest latitudes. 

In the justification of the extraordinary wave ap- 
proximation it was shown that the simplest quasi- 
longitudinal approximation to the phase refractive 
index eq (15) is better than one in common use 
eq (14). 



The work described in this paper was partly sup- 
ported by the International Geophysical Year pro- 
gram of the National Academy of Sciences. 
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